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ABSTRACT 

Much of our understanding of dark matter halos comes from the assumption that the mass-to-light 
ratio (T) of spiral disks is constant. The best way to test this hypothesis is to measure the disk 
surface mass density directly via the kinematics of old disk stars. To this end, we have used planetary 
nebulae (PNe) as test particles and have measured the vertical velocity dispersion (a z ) throughout 
the disks of five nearby, low-inclination spiral galaxies: IC 342, M74 (NGC 628), M83 (NGC 5236), 
M94 (NGC 4736), and M101 (NGC 5457). By using Hi to map galactic rotation and the epicyclic 
approximation to extract a z from the line-of-sight dispersion, we find that, with the lone exception 
of M101, our disks do have a constant T out to ~3 optical scale lengths (Hr). However, once outside 
this radius, <y z stops declining and becomes flat with radius. Possible explanations for this behavior 
include an increase in the disk mass-to-light ratio, an increase in the importance of the thick disk, 
and heating of the thin disk by halo substructure. We also find that the disks of early type spirals 
have higher values of T and are closer to maximal than the disks of later-type spirals, and that the 
unseen inner halos of these systems are better fit by pseudo- isothermal laws than by NFW models. 
Subject headings: galaxies: individual (IC 342, NGC 628, NGC 5236, NGC 4736, NGC 5457) - 
galaxies: kinematics and dynamics — galaxies: spiral — planetary nebulae: general 



1. INTRODUCTION 

Traditionally, spiral galaxies have three main compo- 
nents: a dynamically hot bulge, a dynamically cold disk 
(which can sometimes be divided into "thin" and "thick" 
components), and a mysterious dark matter halo. How- 
ever, measuring the mass of each component separately 
is a challenge, especially in a galaxy's extreme outer re- 
gions where the dark matter halo dominates. Even in 
late-type spirals with minimal bulges, the mass profiles 
of the dark halos cannot be decoupled fro m the visi- 
ble d isk mass using rotation curves alone (|Barnes et al.l 
l2004h . Most analyses have therefore relied on the "max- 
imal disk" hypothesis, wherein the disk mass-to-light ra- 
tio (T) is assumed to be constant with radius, and as 
large as possible to acco unt for the rotation in the inner 
parts of the galaxy (e.g . jKentlll986HPalunas fc Williams! 
&M [Sofu7eTaT]T2003l ). 

Is the disk T really constant and maximal? Perhaps 
the best way to break the disk-halo degeneracy is to mea- 
sure the surface mass of a disk directly from the z mo- 
tions of stars. In systems where the primary stellar mo- 
tion is circular, the phase-space distribution of stars in 
the direction perpendicular t o the plane is fixed by th e 
disk's gravitational potential (jBinney 
For example, in the Milky Way, star counts and veloc- 
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ity measurements toward the Galactic pole initially sug- 
gested tha t the stellar disk is roughly isothermal (e.g., 
IOort][l965h . If this were indeed the case, then the inte- 
grated face-on stellar velocity dispersion, a z , would be 
related to the disk surface mass, S, by 

a 2 z {R) = KGE(R)h z , (1) 

where K — 2tt is a constant, G is the gravita- 
tional constant, and h z is the scale height of the stars 
van der Kruitll l988K Since surveys of edge-on galaxies 



(e.g., Ivan der Kruit fc Searldfj^^TBi^v aev & Mitronova 
2002) imply that h z is roughly constant with radius, mea- 
surements of <j z can produce estimates of disk mass that 
are independent of a galaxy's rotation curve. 

Absorption line studies have shown that, in the in- 
ner regions of spirals, a z generally follo ws the exponen- 
tial law expected from a c onstant T disk (|Bottemalll993l : 
iGerssen et aLlll997l l2000h . However, these surveys were 
limited by surface brightness, and do not extend more 
than ^2 scale lengths from the nucleus. To probe the 
outer regions of disks, where the dark matter halo is 
more important, some other technique is required. Plan- 
etary nebulae (PNe) are ideal tools for this purpose: 
PNe are extremely bright in [O III], abundant out to 
>5 disk scale lengths, measurable to ^2 km s _1 with 
fiber-fed spectrographs, and present in all stellar popu- 
lations with ages between 10 8 and 10 10 yr. Moreover, 
these old disk stars are easy to distinguish from other 
emi ssion-line sources via their distinctive [O III]-Ha ra- 
tio (jCiardullo et al.ll2002fl . 

Planetary nebulae have only recently been used to 
probe the kinematics of spiral galaxies. The largest study 
of this kind has been for the Local Group galaxy M31, 
where the analysis of over 2000 PN velocities has demon- 
strated a flattening of the disk's line -of-sight velocity dis - 
persion at large galactocentric radii (|Merrett et al.l2 006) . 
Although the high-inclination of the galaxy precluded 
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TABLE 1 

Program Galaxies 





Hubble 




Distance 


h R 


Disk 






Number of 


Survey 


Galaxy 


Type a 


i b 


(Mpc) c 


(kpc) d 


Mo e 


E(B — V) ! 


(km S -1 ) s 


PN Velocities 


Region 


IC 342 


Scd 


31° 


3.5 ±0.3 


4.24 


19.60 


0.558 


200 


99 


4.'8 


M74 


Sc 


6.5° 


8.6 ±0.3 


3.17 


20.20 


0.070 


170 


102 


4.'8 


M83 


SBc 


24° 


4.8 ±0.1 


2.45 


19.07 


0.066 


255 


162 


18' 


M94 


Sab 


41° 


4 4+0-1 


1.22 


18.88 


0.018 


200 


127 


5.'8 


M101 


Scd 


17° 


7.3 ±0.5 


4.99 


20.29 


0.009 


250 


60 


8' 



a From RC3 

b IC 342: |Crosthwaite et al.l l|2000f) : M74: IFathi et all 1120071) : M83: IPark et al] ll2001fl : M94: Ide Blok etldl 1120081) : M101: 
IZaritskv et al.l U990D 

c From Paper I, except for M101, which is from Fcldmcicr ct al. (1997); note that this distance has been updated to be 
consistent with M " = -4.47 ICiardullo et al.H200'2T) 

d IC 342: |White &; Bothunl J2003I, J-band); M74: IMollenhofH 1 12004 fl-ba nd); M83: IKuchinski et al.l ( 120001 fl-band); M94: 

Ide Blok et alj l|2008t 3.6/rni) and IGutierrez et al.l d2009t fl-ban d); M101: iKoopmann et al.l (|2006l . fl-band) 

e In mag arcsec " 2 , corrected for galac tic inclination. I C 342: IWhite fc Bothunl (120031 /-band); M74: IMollenhoffl 02004 
fl-band); M83: IKuchinski et al.l d2000l ii-band); M94: IGutierrez et al.l <2009l . fl-band); M101: IKuchinski et al.l (|2000l 
i?-band extrapolation) 

f From[Schlcecl ct U] J1998T 

s IC 342: ICrosthwaite et al.l <2000f) : M74: IFathi et all J2007I) : M83: IPark et al] 11200111 : M94: Ide Blok et"afl 120081) : M101: 
IZaritskv et al.l U990D 



any dynamical measurement of disk mass, these data did 
present evidence for disk flaring at R > 11 kpc. The 
analysis of 140 PN velocities in the moderately inclined 
(i ~ 56°) Local Group spiral M33 yielded a similarly 
surprising result: by assuming a con s tant scale-height 
for the galaxy's disk, ICiardullo et alj (|2004D concluded 
that the disk mass-to-light ratio actually increases with 
radius, going from Ty ~ 0.3 in the inner reg i ons to 
T y ~ 2.0 at ~ 9 kpc. Finally, iDouglas et all pOOOl ) 
measured the velocities of 67 PNe in the low-inclination 
(i ~ 35°) spiral M94, and found a line-of-sight velocity 
dispersion that declines exponentially over ^4 disk scale 

lengt hs . 

Pap ers I and II of this series ([Herrmann "et~aII[2008L 
2009))) presented narrow-band imaging and follow-up 
spectroscopy of PNe in a sample of large (r > 7'), nearby 
(D < 10 Mpc), low-inclination (i < 41°) spiral galaxies. 
Here we use these data to understand the disk kinemat- 
ics of five galaxies: IC 342, M74, M83, M94, and M101. 
We introduce our sample in §2, and in §3, we discuss the 
problem of the disk scale height, the one parameter of 
our analysis that cannot be measured. We explore its 
dependence on galactic radius, its relationship to other 
disk properties, such as scale length, and alternatives to 
the isothermal disk assumption. We then begin our anal- 
ysis in §4 by compensating for the fact that none of our 
nearby galaxies is precisely face-on. We use Hi velocity 
maps to deproject the galactic disks, define the circular 
velocity as a function of radius, and estimate asymmet- 
ric drift. Next, in §5 we identify those PNe which likely 
belong to the galaxies' spheroidal components, and elim- 
inate them from the analysis. We find that very few 
objects fall into this category: only three in IC 342, four 
in M74, six in M83, three in M94, and one in M101. 
Then, in §6, we examine the PNe's line-of-sight disper- 
sions without these halo objects. We show that in the 
inner regions of our galaxies, the line-of-sight velocity 
dispersion, ai os , generally declines as expected for a con- 
stant h z , constant T disk. However, once past ~4 optical 
scale lengths, ai os flattens out at a value much larger than 
expected. We consider a few possible explanations for 
this phenomenon, including the presence of emission-line 



contaminants in our sample, the effects of internal extinc- 
tion, and the behavior of large-scale galactic warps. In 
§7, we model the disks' velocity ellipsoids, and extract a z 
from the line-of-sight velocity dispersions. As expected, 
we find that ur and o$ contribute little to the measured 
dispersion of our nearly face-on disks. In §8, we describe 
the results for each of the five galaxies in our study. We 
convert our values of a z into disk surface masses and cal- 
culate disk mass-to-light ratios. Then, in §9, we analyze 
the results from the inner and outer regions of the galax- 
ies. We compute disk and halo rotation curves, compare 
our disk mass surface densities to those expected from 
maximal disks, show that pseudo- isothermal cores fit the 
residual rotation curves better than NFW models, and 
discuss the physical mechanisms that may produce a con- 
stant a z , including the interaction of thin and thick disks 
with halo substructure. We conclude in §10 by summa- 
rizing our results and suggesting future observations that 
will assist in their interpretation. 

2. THE SAMPLE 

Table 1 lists the five galaxies observed in this sur- 
vey, along with their physical properties. IC 342, M74, 
and M101 are large, late-type spirals with little or no 
spheroidal components, and rotation velocities greater 
than 170 km s -1 . M94 is an earlier system (type Sab) 
with a small, but bright bulge; in this galaxy, our PN 
identifications were largely confined to regions outside of 
~1.5 kpc, where the disk of the galaxy dominates. Fi- 
nally, M83 is a strongly barred SBc system with isophotal 
twists and a warped outer disk. 

Our galaxies were initially surveyed for planetary neb- 
ulae with on-band/off-band [O III] and Ha imaging with 
the WIYN telescope in the north and the CTIO Blanco 
telescope in the south (Paper I). Follow-up spectroscopy 
was then accomplished via the Hydra multi-object spec- 
trographs of these same two telescopes (Paper II). These 
observations produced between 60 and 162 PN radial ve- 
locities per galaxy, with a median measurement error of 
a v ~ 6 km s , and, in all cases, a v < 15 km s . The 
positions, magnitudes, radial velocities, and velocity un- 
certainties of all the PN candidates are given in Paper II. 
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TABLE 2 
Scale heights (in pc) 



Galaxy 


dG a 


K02 b 


BM02 C 


Initial Range 


Stability d 


Rotation Curve e 


Final Range 


IC 342 


505 


424 


565 


400 - 600 


> 300 


> 100 


400 - 600 


M74 


369 


327 


423 


300 - 500 


> 300 


> 170 


300 - 500 


M83 


285 


377 


402 


200 - 500 


> 300 


> 200 


300 - 500 


M94 


277 


186 


218 


100 - 300 


> 300 


> 200 


300 - 400 


M101 


623 


523 


697 


500 - 700 


> 300 


> 60 


500 - 700 



a Estimated using hn/h z values from dc Griis ( 19 981) 
Estimated using hn/h z values from|Krcgcl ct al. (200^) 

c Estimated using hn/h z values from Bizyacv & Mitronova (2002) 
Stability limit on h z 

c Rotation curve limit on h z 



3. ESTIMATING SCALE HEIGHT 

Before we can translate our PN velocities into disk 
mass, we need to adopt some form for the phase-space 
distribution of stars perpendicular to the disk. This is 
not an entirely settled question, even in the Milky Way. 
As described by Ivan der Kruitl |l988), the vertical struc- 
ture of stars in a galactic disk can be parameterized by 
a family of models 



p(z) oc sech 2/ ™ ( — 



(2) 



which have isothermal (n = 1) and exponential (n = oo) 
distributions as their limiting cases. Depending on the 
value of n, Poisson's equation yields different constants 
for the relationship between the disk mass surface density 
and the integrated face-on velocity dispersion. In other 
words, in equation (Q]), 

n = 1 => Piso(z) oc scch 2 (z/2/i 2 ) => K = 2n (3) 
n = 2 => Pim(z) oc sech(z/h z ) => K = 1.70517T 
n = oo Pcxp(z) oc exp(— z/h z ) => K = 1.5ir. 

IVan der Kruitl (|1988l ) has argued that the intermediate 
(pi m ) case best fits star counts in the Milky Way, espe- 
cially at low z, where the superposition of many stellar 
populations, each with their own value of a z , blend to- 



;ether . To take this blending into account. ISiebert et al.l 
2003) combined several isothermal components to find 
a result similar to the exponential case. However, when 
initial m ass functions an d star formation rates were con- 
sidered, lHavwood et alj (|1997l ) found that a single ex- 
ponential was insufficient to model the Galaxy's verti- 
cal st ructure. More recently, [Schwarzk opf fc Dettmarl 
( 2000) have shown that in a sample of 61 edge-on spi- 
rals, 60% have vertical distributions that are best fit by 
the sech(z) law. However, unlike the Milky Way star 
counts, these measurements were the result of a fit to sur- 
face brightness, and therefore represent a solution that 
is weighted by luminosity (actually, by the inverse of the 
stellar mass-to-light ratio), rather than the actual mass 
density of stars. 

While PN observations do not strictly sample mass, 
neither do they produce results that are luminosity- 
weighted. Generally speaking, PNe are excellent trac- 
ers of light, but only in populations where the mass-to- 
light ratio chang es slowly with a ge, i.e., in systems older 
than ~ 0.1 Gyr (|Maraston| [l998'l . Younger populations, 
which may dominate the luminosity of a galaxy, produce 
no planetary nebulae whatsoever. Thus, the distribution 



of PNe should be closer to one that is space-weighted, 
rather than one that is luminosity weighted. 

Given these issues we chose to model our disks using an 
intermediate (n — 2) vertical distribution. This assump- 
tion increases our derived disk masses by 17% over the 
isothermal case, while decreasing it by ~ 14% compared 
to exponential distributions. 

With the vertical distribution set, we next needed to 
establish the values of h z appropriate for our sample 
of galaxies. Since all of our targets are low-inclination 
systems, their vertical scale heights are not directly ob- 
servable. Consequently, in order to constrain h z , we 
had to use correlations with other measurable quantities, 
such as disk scale length (ha) and disk surface bright- 
ness. Several studies have investigated these trends, 
generally finding h R /h z ~ 10 for large, non-interacting 
late- type spirals, such as those observed in our sur- 
vey ( e.g.. Ide Griis fc van der Kruitl 119961 : iKregel et al. 
2002t lYoachim fc Dalcantonll2006lh ~ In addition, Ide Griisl 
(|1998h and lKregel et al.l (|2002f ) have explored the behav- 
ior of h]i/h z versus Hubble type. Both studies find that, 
though later-type spirals generally have larger h z values 
than earlier type objects, they also have larger values 
of hji/h z . Thus, while Scd galaxies have hn/h z ~ 10, 
the Sab systems have rati os as low as ~5. Finally, 
iBizvaev fc Mitronova! (|2002f ) have examined the depen- 
dence of hn/hz as a function of -ftf-band central surface 
brightness, K , finding 



h R /h z = (0.65 - 0.0295ifoV 



(4) 



This relation has a substantial amount of scatter, but if 
we use the radial surface brightness profiles of the 2 Mi- 
cron All Sky Survey (2MASS) Large Galaxy Atlas web- 
site 4 to estimate Kq, then the implied values of hn/h z 
for the galaxies in our sample range from ^5.6 to ~7.5. 

Table [2] summarizes our initial estimates for h z , based 
on the galactic scale lengths given in Table 1. In §7, we 
will further investigate the plausibility of these numbers 
using stability arguments applied to the observed PN 
kinematics. We will also consider whether the use of 
a single value of h z is appropriate, or whether "flaring" 
must be considered in the extreme outer disks of galaxies. 

4. VELOCITY FIELDS FROM PLANETARY NEBULAE 

The line-of-sight velocity (vi os ) of a planetary nebula 
rotating in the disk of a galaxy is given by 

vios = cos (j) sin i + Vr sin (f) sin i + v z cosi + v sysi (5) 

4 http: / / www.ipac. caltech. edu / 2mass / gallery /largegal /atlas . html 
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FlG. 1. — [O III] A5007 images of our program galaxies, with the positions of the PNe superposed. Left: observed PN radial velocities. Middle: 
VLA Hi velocities (outer annuli) and our model of galactic rotation (inner dots). Right: velocity residuals formed by subtracting the VLA Hi 
velocities (corrected for asymmetric drift) from our PN velocities. Color coding for all but the M74 Hi figure is in 30 km s — 1 increments from 
> 60 km s -1 systemic (red) to < —60 km s _1 systemic (violet). Hi velocities for the face-on M74 system are shown in 10 km s _1 increments, 
from 20 < v < 30 km s — 1 to —30 < v < —20 km s — . Solid points are well-determined velocities ((7 V < 15 km s — 1 ), crosses are uncertain 
velocities {a v > 15 km s — 1 ), squares represent probable halo objects, and small wh ite dots denote photo metrically identified PNe with no velocity 
measurement. The green boxes in the left panel of M94 show the fields studied bv lDouglas et al.l J2000D . The axes in the middle panels illustrate 
our adopted values for the galactic position angles. North is up and east is to the left. 
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where v sys is the galactic systemic velocity, i the disk in- 
clination, and vr, Vij), and v z are the PN's component 
velocities in the radial (i?), azimuthal (cf>), and verti- 
cal (z) directions. Note that even though our galaxies 
are nearly face-on, the effect of rotation is still consider- 
able. This is easily seen in the left panels of Figure [TJ 
which display our measured PN velocities superimposed 
on an [O III] A5007 image of each galaxy. Clearly v$ 
is much larger than v z even when it is diminished by 
sini. (The one exception is M74 which has an incli- 
nation of only ~6°.) Figure [2] shows the rotation in 
another way by plotting our measured velocities as a 
function of azimuthal angle, 9 (north through east in 
the plane of the sky) . The location of the maximum of 
each sinusoid indicates the position angle of the princi- 
pal axis, while the amplitude of the curve places a limit 
on v max smi. This rotation must be removed before o~ z 
can be isolated. Fortunately, The H i Nearby Galaxy Sur- 
vey (THINGS: IWalter et al.ll200ll collaboration has re- 
cently used the Very Large Array (VLA) to create very 
high resol ution velocity map s of fo ur of our five spirals. 
Moreover. I Crosthwaite et alj (|2001f ) have made a similar 
Hi map for IC 342, though with lower resolution. 

4.1. Removing Galactic Rotation 

To remove galactic rotation from our PN velocities, we 
first needed to define each disk's three dimensional ge- 
ometry, and deproject the coordinates from the plane of 




50 100 150 200 250 300 350 
8 (Angle East of North in the sky, in degrees) 



FlG. 2. — PN radial velocities as a function of position angle on 
the sky, where 9 is the angle East of North. The sinusoidal pattern 
illustrates galactic rotation, and the location of maximum velocity, 
Vmax sinz, defines the position angle of the major axis. Solid points 
are well-determined velocities (cr v < 15 km s — 1 ), crosses are uncer- 
tain velocities (a v > 15 km s _1 ), and squares represent probable halo 
objects. 



TABLE 3 

Geometric Parameters for Disk Models 



Galaxy 


PA 


i 


Radial Range 




(in degrees) 


(in degrees) 


(in arcmin) 


IC 342 


45 


31 


Full range 


M74 


25 


9 


Full range 


M83 


226 


21 


R < 5.8 




246 


27 


5.8-8.0 




236 - 226 


29 - 30 


8.0 - 9.0 a 




226 - 201 


30 - 35 


9.0 - 14.0 a 




201 - 172 


35 - 46 


14.0 - 25.0 a 


M94 


305 


37 


R < 2.45 




295 


37 


R > 2.45 


M101 


35 


17 


Full range 



a PA and i are linearly dependent on R in these regimes. 
(Based on model in jHuchtmeier fc Bohn cnstcngcl ( 19811) .') 



the sky to the plane of the galaxy. We started with liter- 
ature values for each galaxy's position angle, inclination, 
and rotation curve, and then iteratively refined the pa- 
rameters to fit the VLA Hi velocity maps. In some cases, 
this meant changing the quoted position angle of the ma- 
jor axis by a few degrees; at other times, the inclination 
or rotation amplitude required modification. When the 
latter condition occurred, we preferred to tweak the in- 
clination by less than five degrees rather than adjust the 
rotation velocity by 20 to 60 km s _1 . Table [3] presents 
our final geometric parameters for our disks. 

The middle panels of Figure [1] compare our final mod- 
els for galactic rotation with the measured Hi velocity 
at the location of the objects tabulated in Paper II. For 
IC 342, M74, and M101, a simple thin disk with a con- 
stant position angle and inclination works well over the 
entire region of our study. For the more extended survey 
of M94, two different position angles were needed, one 
for the inner disk (PA = 305°), and one for the system's 
outer ring (PA = 295°). Finally, to model the known 
isophota l twists of M83, we started with the w arped disk 
model of lHuchtmeier fc Bohncn stengell (|1981l ). which ex- 
tends well beyond our survey area, and then adjusted the 
parameters to fit the higher resolution THINGS data. 
The axes in the middle panels of Figure [1] illustrate our 
final values for the galactic position angles. 

We do note that just because our Hi residuals exhibit 
no systematic trend with azimuth or radius, that does 
not mean that a map of the stellar velocity would be 
just as wel l-behaved. However, it is likely that the two 
are linked. iMorrison et al.l (| 19971 . and references therein) 
have found that disks which contain warps in Hi often 
have stellar warps as well. Moreover, in at least one 
case (NGC 5907), bo th warps have similar orientations 
(jMorrison et al.lll994D . Thus, our procedure of modeling 
the disk orientation using Hi data should be valid. 

4.2. Estimating Asymmetric Drift 

Once the geometry of the disk was set, we subtracted 
off the Hi rotation from the observed PN velocities. This 
generated the velocity residuals shown in Figure [3l As 
the figure illustrates, the Hi rotation curves do not re- 
move all trace of systematic behavior from the veloci- 
ties. In particular, the PN velocities of M83 and M94 
now show a low-amplitude convex bowing with respect 
to 180° azimuth. This is the signature of asymmetric 
drift, i.e., the lag between the circular velocity, as mea- 
sured by the gas, and the average azimuthal velocity of 
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TABLE 4 

Asymmetric Drift Estimates (km s _1 ) 



50 100 150 200 250 300 350 
(Azimuthal angle in the galaxy plane, in degrees) 

FlG. 3. — PN velocities after the removal of disk rotation, as defined 
by the motion of the Hi gas. Note that for MS3 and M94, the pat- 
tern of velocity residuals has a slightly convex appearance. This is the 
signature of velocity oversubtraction, due to the action of asymmetric 
drift. Sec Figurc[2]for an explanation of the point types. 

the stars. By using the Hi data to define rotation, we 
oversubtracted the stars' mean azimuthal motion. 

To estimate the amplitude of this effect, we can use the 
velocities of the PNe themselves. Because galactic disks 
are cold systems, the i^cos^sini term of equation ^ 
clearly dominates. Thus, to a first approximation, 



Vic 



J sys 



C( )S Sin I 



(6) 



Obviously, at tfi = 90° and cf> = 270° the solution to this 
equation is undefined, but if we exclude PNe within 10° 
of the minor axis and bin by radius, we can derive mean 
values of v$ over the entire disk of the galaxy. These 
values can then be compared to the observed motion of 
the gas. This analysis confirms the results of Figure [3J 
that asymmetric drift is negligible in IC 342 and M101, 
~15 km s" 1 in M74, between 20 and 40 km s" 1 in M83, 
and between 40 and 50 km s _1 in M94. 

Are these velocities reasonable? It is important to note 
that asymmetric drift is extremely difficult to measure, 
so very few estimates exist in the literature. In the solar 
neighborhood, where the rotational velocity is v max ~ 
220 km s -1 , the asymm e tric drift is v as d ~ 20 km s _1 
(jDehnen fc BinnevNl998t lOlling fc Penned 12003ft. This 
is rou ghly twice the value found bv ICiardullo et alJ 
(2004) from radial velocity measurements of PNe in 
the small, late- type spiral M33 (v max ~ 100 km s _1 ). 
For the large Sb spiral M31 (v max ~ 260 km s _1 ), 
PN measu rements yield value s of v asl j between 40 to 
50 km s _1 (M errett et al.|[2006T ). and in the barred lentic- 
ular NGC 1023 (v max ~ 220 km s" 1 ), v asd ~ 40 km s" 1 
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Estimate 


^max 




0. 2f ma:E 


Best Value 


IC 342 


negligible 


200 


20 


40 


20 


M74 


15 


170 


17 


34 


34 


M83 


20 - 40 


255 


26 


51 


51 


M94 


40 - 50 


200 


20 


40 


40 


M101 


negligible 


250 


25 


50 


25 



(jNoordermeer et "al"1l2008l ). These estimates suggest that 
asymmetric drift is usually 10-20% of the rotation ve- 
locity in late-type galaxies and slightly larger in earlier 
systems. Both correlations make sense: in the Milky 
Way, asymmetric drift is correlated with stellar popu- 
lation, since older stars have had more time to scatter 
(jDehnen fc Binnevlll998| ). and the amount (and ampli- 
tude) of this scatter depends on the speed of rotation. 

Table |4] compares our observed asymmetric drift ve- 
locities to values which might be inferred from a simple 
scaling of the rotation curve. To facilitate this compar- 
ison, we have used a single drift value for each galaxy; 
this is consistent with the PN data of M31, M33, and 
NGC 1023, which suggest that the asymmetric drift 
within a galaxy is roughly constant with radius. Wc 
also list the "best fit value" for each object, i.e., the one 
that minimizes our residual line-of-sight PN velocity dis- 
persions while still being in the range between 0.1 and 
0.2v max (see next section). As Figure[3]and Table2]indi- 
cate, our galaxies exhibit a range of drift velocities. M83, 
which is a strongly barred galaxy and the fastest rota- 
tor in our sample, has the largest value of v as d- In M94, 
which is an early- type system (Hubble type Sab), v as d is 
also rather large, ^40 km s _1 . On the other hand, IC 342 
and M101 both have rather low values for the asymmetric 
drift. This suggests that their disks are relatively young, 
an idea supported by the fact that these two galaxies 
are the latest systems in our sample (Hubble types Scd). 
Finally, the asymmetric drift for M74 is slightly higher 
than expected; however, due to the galaxy's extremely 
low inclination (~ 6°), this number is very uncertain. 
We will revisit our estimates of v as d in §9, when we have 
measurements for the disks' radial and perpendicular ve- 
locity dispersions. 

The right panels of Figure Q] display our PN line-of- 
sight velocities with the Hi gas rotation and our best- 
fit asymmetric drift values removed. The data show no 
trend with cj>, indicating that the adopted drift values are 
reasonable. We note that the results that follow depend 
very little on the exact values of v as d- because our galax- 
ies all have low-inclinations, 10 to 20 km s _1 changes in 
the asymmetric drift have almost no effect on the final 
values for the z component of the velocity dispersion. 

5. IDENTIFICATION OF HALO OBJECTS 

After the removal of rotation, our next step was to 
exclude those objects not part of the general disk pop- 
ulation. Since all five of our program galaxies are disk- 
dominated systems, the vast majority of the PNe should 
be rotating close to the galactic plane. However, a 
few halo objects may be present, and these objects can 
greatly distort our measurement of the disk velocity dis- 
persion. For example, a halo PN will typically orbit at 
a velocity close to v max , but with some random orien- 
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IC 342 Galactic Radius (R-band scalelengths) 

FlG. 4. — IC 342 PN radial velocities (top), linc-of-sight velocity 
dispersion (middle), and velocity dispersion in the z direction (bottom), 
all as a function of galactic radius. PNe rejected as non-disk objects are 
identified in the top panel, along with the amount of their departure 
from the mean disk; these objects are not included in the lower plots. 
PNe with velocity uncertainties greater than 15 km s — 1 are not shown; 
some Hll region contaminants arc plotted as light gray crosses. The 
light gray triangles and black curve in the top panel indicate the average 
residual velocity in each bin. The alternating black and dark gray 
points illustrate the radial bins. In the two lower panels, the dotted 
line shows our median velocity error, the solid black curve denotes an 
exponential with twice the galaxy's /-band scale length, and the gray 
curve adds these two relationships in quadrature. This last curve is 
the one expected for a constant h zl constant T disk. The radial width 
of each 16 object bin is illustrated via the errors bars in x\ \a velocity 
uncertainties arc indicated by the errors in y. For IC 342, the velocity 
dispersion agrees with the simple model of a constant h z , constant T 
disk. 

tation in the galaxy. For a nominal rotation speed of 
^200 km s~ x , this implies a velocity of ^100 km s _1 
with respect to the local standard of rest. This is much 
larger than the expected value for velocity dispersion of 
disk stars. 

Still, identifying halo PNe can be difficult. Depend- 
ing on the orientation of its orbit, a halo PN can have a 
line-of-sight velocity that mimics that of the disk. More 
importantly, if the disks of our program galaxies are any- 
thing like that of the Milky Way, we can expect a line-of- 
sight velocity dispersion of ~40 km s _1 . Thus a typical 
halo planetary may have a velocity that is only ~2.5 a 
away from the mean of the disk population. 

Nevertheless, we can identify obvious halo candidates 
via a sigma-clipping algorithm. For each PN, we iden- 
tified its 15 closest neighbors (in radial distance) and, 
using those objects, we estimated the local disk velocity 
dispersion. Any PN more than 2.5 standard deviations 
away from the mean (where a is defined from the disk 



dispersion and PN velocity uncertainty added in quadra- 
ture), was flagged as a possible halo contaminant and 
eliminated from the analysis. Formally, this means that 
out of our entire sample of ~550 PNe, we may be im- 
properly excluding as many as ~11 disk objects. How- 
ever, given the deleterious effect that halo PNe may have 
on our results, we feel that 2.5 a is a reasonable threshold 
for PN exclusion. 

The top panels of Figures [4][8] show the velocity resid- 
uals of all our PNe and identify those objects with dis- 
crepant velocities. Even with our harsh 2.5 a condition, 
very few PNe are excluded from the analysis. In M74, 
PN 127 is an obvious (-6 a) outlier and PNe 66, 96, 
and 146 also appear to be contaminating halo objects. 
In IC 342, PNe 42, 98, and 111 are likely halo interlop- 
ers, while in M83 PNe 8, 43, 48, 50, 102 and 175 have 
the kinematics of the spheroidal component. Despite the 
fact that M94 is an Sab galaxy, only three of its PNe 
(26, 63, and 102) are questionable: this is likely due to 
the fact that the galaxy's bulge is small, and that our 
PN identifications were largely confined to the galaxy's 
outer regions where the disk dominates. Finally, none 
of the PNe in MIOl's kinematic sample have discrepant 
velocities: although PN 51 is a possible halo object, we 
had already excluded it from the analysis due to its large 
measurement error. 

M74 Galactic Radius (kpc) 
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M74 Galactic Radius (R-band scalelengths) 

FlG. 5. — M74 PN radial velocities (top), line-of-sight velocity disper- 
sion (middle), and velocity dispersion in the z direction (bottom), all 
as a function of galactic radius. See the caption of Figurc[^for more in- 
formation. Each radial bin contains 16-17 objects. As with IC 342, the 
dispersion curve agrees with that expected for a constant h z , constant 
T disk. Our value of 53 ± 5 km s _1 for the central z velocity dispersion 
is consistent with the value of ~60 km s -1 der ived from integral-fi eld 
spectroscopy of the system's central kiloparscc ( Ganda ct al. 2006). 
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M83 Galactic Radius (R— band scalelengths) 

FlG. 6. — M83 PN radial velocities (top), linc-of-sight velocity dis- 
persion [middle), and velocity dispersion in the z direction [bottom), 
all as a function of galactic radius. See the caption of Figure ^ for 
more information. Each radial bin contains 15-16 objects. The gray 
line in the lower panel illustrates our best-fit two-component model for 
the galaxy. M83 has a strong bar, and the small values of cr; os and 
a z in the second and third bins may be due to the bar's outer Lind- 
blad resonance, which is located at Ft. ~ 7 kpc. If these two points are 
excluded from the analysis, then the dispersion profile looks similar 
to that of M94, with an exponentially declining inner region, and an 
asymptotically fiat outer disk. 

6. RESIDUAL VELOCITIES AND THE LINE-OF-SIGHT 
DISPERSION 

With rotation removed and the halo PNe eliminated 
from the analysis, we could examine the kinematics of 
the galactic disks. We started by assigning an equal 
number of PNe per radial bin (16 per bin in IC 342, 
M74, and M83, 18 per bin in M94, and 15 per bin in 
M101), and computing the line-of-sight velocity disper- 
sion within each bin. According to equation ([1]), this 
function should be a declining exponential with a scale 
length twice that of the light if the following conditions 
hold true: 

1 . The uncertainties in our velocity measurements are 
small compared to a z . 

2. The line-of-sight velocity residuals are dominated 
by the motions of stars in the z direction. 

3. The scale height, h z , is constant over the region of 
the survey. 

4. The disk mass-to-light ratio, T, is constant over 
the region of the survey. 

5. The scale length, h,R, is constant over the region of 
the survey. 



The middle panels of Figures [4][8] compare our disper- 
sions to such an exponential. It is immediately obvious 
that in IC 342 and over the inner ~3 scale lengths of 
M74 and M94, the measured velocity dispersion falls off 
exactly as expected for a constant T disk. This result 
is consistent with that found from integrated-light ab- 
sorption line spectroscopy, although the latter technique 
has never prob ed radii larger than ~2 scale lengths (e.g., 
lBottemal [l993). However, in the outer bins of M94, the 
velocity dispersion does not fall off. Instead it asymp- 
totes out at ^30 km s _1 , a value that is an order of mag- 
nitude greater than the ~3 km s _1 median error of the 
data. This behavior is also present in the outer regions 
of M83: if one ignores the low dispersions of the second 
and third bins (which are affected by the galaxy's outer 
Lindblad resonance), then the dispersion profile is con- 
sistent with a declining exponential in the inner regions, 
and a constant value at large radii. Only M101 does not 
follow the pattern: its line-of-sight velocity dispersion is 
almost flat over ~4 optical scale lengths. 

Before discussing this behavior further, we consider 
four possible systematic errors which could be associated 
with our measurements. 
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M94 Galactic Radius (R-band scalelengths) 

FlG. 7. — M94 PN radial velocities (top), line-of-sight velocity dis- 
persion (middle), and velocity dispersion in the z direction (bottom), 
all as a function of galactic radius. See the caption of Figure^for more 
information. Each radial bin contains 17-18 objects. The gray line in 
the lower panel illustrates our best-fit two-component model for the 
galaxy, using an internal scale height of 300 pc. Although the system's 
Z velocity dispersion initially falls off with the light, the profile flattens 
dramatically in the outer disk. 
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M101 Galactic Radius (R-band scalelengths) 

FlG. 8. — M101 PN radial velocities (top), line-of-sight velocity 
dispersion (middle), and velocity dispersion in the z direction (bottom), 
all as a function of galactic radius. See the caption of Figure ^ for 
more information. Each radial bin contains 15 objects. The gray line 
in the lower panel shows our best-fit exponential. Note that MIOl's 
disk velocity dispersion is roughly flat with radius. 

6.1. Observational Uncertainties 

Measurement errors provide a lower limit to our ob- 
served disk velocity dispersion. If these errors were large, 
it could produce a plateau in the observed dispersion 
curve. However, as demonstrated in Paper II, the uncer- 
tainties quoted for our PN velocities are accurate, and 
significantly smaller than the measured disk dispersion 
(see the dotted lines in Figures HUH]). Alternately, if some 
intrinsic property of the PNe themselves limited the pre- 
cision of our measurements, we would again see a plateau 
in o~ z . For example, if the measured radial velocity of an 
(asymmetrical) PN depended on its orientation, then we 
would expect the dispersion profile to have an artificial 
floor. Of course, any such floor would need to be less 
than the ~12 km s _1 value recorded in the outer disk of 
M74. Moreover, it is difficult to conceive of how a for- 
bidden line with a ^20 km s" 1 symmetric velocity width 
could yield a ~12 km s _1 velocity error which depends 
on viewing angle. In any case, since this floor should be 
the same in every galaxy, it cannot explain dispersions 
of 20 to 30 km s" 1 seen in M83, M94, and M101. 

One other possible error comes from the small 
(~5 km s _1 ) uncertainties associated with the Hi mea- 
surements used to remove rotation, and with our esti- 
mated values for asymmetric drift. Such an error is sup- 
pressed by sini, but it can still affect our values for ct; os . 
To investigate this problem, we performed a series of ex- 
periments, in which the asymmetric drift was allowed to 



vary between and 50 km s _1 . These simulations showed 
that our line-of-sight dispersion measurements are insen- 
sitive to errors in the global rotation, changing by only a 
couple of km s _1 . This type of error cannot be responsi- 
ble for the dispersion profiles observed in the figures. 

6.2. Emission Line Contaminants 

Papers I and II describe the criteria used to identify the 
planetary nebulae used in this analysis. The data were 
well scrutinized, and any bias in the selection process 
was made on the side of eliminating valid PNe rather 
than including possible contaminants. Nevertheless, it 
is still conceivable that compact, high-excitation Hil re- 
gions have slipped into the kinematic sample. The most 
serious case is for M101, where we have included six 
objects whose spectra may contain some contamination 
from a kinematically distinct source, and six other ob- 
jects with questionably bright Ha lines. However, even 
without these objects, MIOl's line-of-sight velocity pro- 
file is flat with radius. 

In addition, during our M74 and M83 observations, we 
obtained spectra of a significant number of suspected Hil 
regions. To examine the effect these objects can have on 
our results, we repeated the analysis of these two spirals, 
while including 21 and 63 contaminating sources, respec- 
tively. As expected, the addition of a younger population 
decreases the observed dispersion. In M74, the new, con- 
taminated oios was smaller, but within the uncertainties, 
still consistent with the results from the uncontaminated 
data. In M83 where we added three times the number of 
contaminants, the dispersion again decreased, especially 
in the inner two bins and the outermost bin, where most 
of the interlopers are located. Nevertheless, the general 
shape of the dispersion profile remained the same. The 
results confirm that the addition of a young Population I 
component can only decrease the dispersion; it cannot 
explain the relatively high values of o~i os in the outer 
disks of our galaxies. 

6.3. Internal Extinction and Selection Effects 

Since spiral galaxies are known to have dust at low 
galactic latitudes, one might wonder if internal extinc- 
tion could affect the results of our survey. Practically 
speaking, a thick layer of dust would preferentially dim 
those PNe near the galactic plane and on the far side of 
the galaxy, while leaving objects with large z distances 
on the near side unaffected. In the case of an isothermal 
disk, a z is independent of z, so dust would have no ef- 
fect on our results. However, for intermediate and expo- 
nential disks, the observed v elocity dispersion in creases 
with height above the plane (jvan der Kruitl[l988f ). Con- 
sequently, a thick layer of dust could theoretically drive 
dispersion measurements towards higher values. How- 
ever, such a model is unlikely, as it is inconsistent with 
the number of PNe we observe (since the dust would ex- 
tinguish large numbers of PNe below our detection limit), 
and the distribution of PN [O III] /Ha line ratios (see Pa- 
pers I and II). Moreover, because the radial scal e length 
of du st is usually larger than that for the stars (|Bianchil 
l2007f) . it is the dispersion me asurements in the inne r disk 
that would be most effected. ICiardullo et al.l (|2004l ) have 
discussed the situation in more detail and found that 
dust cannot explain the flat velocity dispersions at large 
galactic radii. 
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6.4. Incorrect Disk Models and Large Scale Warps 

Warps are known to exist in the outer regions of spiral 
galaxies and an undetected change in inclination could 
propagate into incorrect values of oi os . In fact, one of the 
main reasons we adjusted the geometric parameters when 
fitting the VLA Hi velocity maps was to remove warping 
as a cause for concern. As indicated in §4.1 and the 
central panels of Figure [U only M83 and M94 required 
anything other than a simple, flat disk. However, since 
these are the best examples of flat dispersion profiles, it 
is important to re-examine this possibility. 

M94's fit required two components, one for the galaxy's 
interior regions (R < 5 kpc), and one for the outer disk, 
so it is unlikely that our warping model could be an issue 
in this system. The model for M83, however, was based 
on pre viou s Hi analyse s bvlHuchtmeier fc Bohnenstengell 
((l98lh and lPark et all (120011) . and was much more com- 
plicated. Yet despite this complexity, warping only 
slightly alters the dispersion measurements and the over- 
all effect on o z (R) is minimal. Specifically, if we assume 
a simple, flat disk, and re-run the analysis of M83's PNe, 
the results are essentially the same. Even a moderate 
amount of warping cannot be responsible for the veloc- 
ity dispersions seen in the figures. 

7. MODELING THE VELOCITY ELLIPSOID 

The velocity dispersions shown in the middle panels of 
Figures[5]through[8]suggest that the outer disks of galax- 
ies may not be simple constant scale height, constant T 
exponentials. However, to confirm this result, we would 
like to extract o~ z from the line-of-sight measurements. 
To do this, we began by breaking down the observed ve- 
locity dispersion into its three orthogonal components: 
the radial (or), azimuthal (er^) and perpendicular (o z ) 
velocity dispersions, i.e., 



'los 



o\ cos 2 . 



> sin 2 i - 



o R sin 



> sin 2 i + a 2 cos 2 i + cr 2 



(7) 

(Note that the (known) observational uncertainty of each 
velocity is explicitly included via the <x meas term.) We 
then removed o^ from the equation by taking advantage 
of the fact that the motion of stars in a cold disk can 
be approximated by the superposition of a circular orbit 
with movement in a small ellipse with frequency 



k =r( 2 + 2 ^r. 



1/2 



(8) 



where V c is th e circular velocity inferred from the motions 
of the Hi gas ()Binnev fc Tremainelll987t) . The azimuthal 
and radial velocity dispersions are therefore related by 
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The line-of-sight dispersion can then be written as a func- 
tion of two variables: or, which modulates with position 
angle in the galaxy, and o z , which does not. We can 
therefore determine the likelihood that a particular com- 
bination of these two variables will generate an observed 
set of PN velocities. 

Figures [9lfT31 display the relative probabilities of such 
an analysis. For the calculation, we have constrained 
the ratio of the two dispersions to lie in the range 
0.25 < o z /or < 1-0, where the upper limit comes 



from the physics of disk scattering, sinc e both theory 
(lVillumsenlll985t Uenkins fc Binnevlll990l) and observa- 
tions ( Bienavmelll999HGerssen et al.l fT997l[2000h require 
or to dominate. The lower limit is a consequence of the 
buckling or "fire hose" instability, which can occur in a 
cold disk (|Toomrdll96l : lAraklfl98a iMerritt fc Sellwoodl 
1994). These limits are responsible for the cone-like 
shapes of the diagrams. 

As the figures show, the shapes of the likelihood prob- 
abilities vary from galaxy to galaxy. This reflects differ- 
ences in galactic inclination. In M74, which is extremely 
close to face-on (i ~ 6°), oi os is almost entirely due to o z , 
and or is essentially unconstrained. Conversely, as the 
disks become more inclined, our sensitivity to ctr grows 
at the expense of our measurement of o z . 

The bottom panels of Figures marginalize these 
likelihoods over radial dispersion to find the most likely 
values of o z and their 1 o uncertainties. Once again, 
the black curves indicate exponential laws with twice the 
scale length of the i?-band light. As expected for low- 
inclination systems, the o z profiles track those of the 
line-of-sight dispersions extremely well. The ratio of o z 
to oi os is roughly constant across the face of each galaxy, 
ranging from -0.76 in IC 342 to -0.86 in M74. For 
IC 342 and M74, the disks are well-fit by a constant 
scale height, constant T model. In contrast, the data for 
M83 and M101 do not follow the expected exponential: 
in these objects, either the scale length, scale height, or 
mass-to-light ratio is changing. Finally, M94 appears to 
be a hybrid: while the inner ~4 scale lengths can be 
fit with a simple model, the outer regions demand an 
additional term. 

One additional constraint can be applied to the data. 
In order for our thin stellar disks to be stable ag ainst ax- 
isymm etric perturbations, they must satisfy the iToomrel 
(1964J) criterion 

3.36GE 

o R > . (10) 

K 

If we combine this relation with equation {1} , we obtain 
a new constraint on the possible values of or and o z 

/Kh zK QR\^ 

V 3.36 J v ' 

where K — 1.7051-/T for intermediate disks, h z is the stel- 
lar scale height, and k is the epicyclic frequency given by 
equation ©. 

The solid curves in Figures [91 1131 show the maximum 
values of o z allowed by the Toomre criterion, given sev- 
eral different estimates of h z (see Table [5]). For example, 
in IC 342, the assumed range of h z is consistent with the 
derived likelihood probabilities; if we were to hypothesize 
a much lower value for h z , then the likelihood of a stable 
solution would decrease dramatically. Similarly, for most 
of M74 and M101, the probabilities lie safely below the 
Toomre curves in the region of stability. However, in the 
galaxies' outermost regions, slightly larger scale heights 
may be in order. Finally, the probabilities for M83 and 
M94 suggest that our assumed values of h z are much too 
small, especially at large galactic radii. This result may 
not be too surprising, since these are also the objects 
that have the most noticeable asymmetric drift. It is 
likely that the energy missing from disk rotation has, at 
least in part, been translated into larger z motions for 
the stars, and a greater galactic scale height. 
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TABLE 5 
Disk Mass Models 



Galaxy 


Filter 


A»(0) 


h z (pc) 


h R (kpc) 


a z (0) (km s" 1 ) 


E(0)(M Q pc' 2 ) 


T(0) 


dlog(Th z )/dr 


IC 342 


/ 


19.60 ±0.20 


500 


4.24 


41 ±6 


146 ± 45 


0.25 ±0.10 





M74 


R 


20.20 ±0.05 


400 


3.17 


53 ± 5 




1.4 ±0.3 





M83 


R 


19.07 ±0.05 


400 


4.0 


73 


578 


1.0 ±0.5 











1200 


10 


40 






0.16 


M94 


R 


18.88 ±0.05 


300 


1.22 


91 ± 15 


1200+f 7 ° 


1.8 ±0.4 







R 


22.69 ±0.10 


900 


13.2 


29 ± 5 


41 ± 15 


2.0 ±0.4 





M101 


R 


20.29 ±0.10 


600 


24.8 


38±| 


ioo±i° 


0.6 ±0.3 


0.07 



8. INDIVIDUAL GALAXIES 

Table [5] summarizes our results for the five galaxies in 
our program. Each object is discussed below. 

8.1. IC 342 

IC 342 is a large (i? > 20 kpc), nearby (D ~ 3.5 Mpc) 
Scd spiral that overfills the ~10' x 10' field-of-view of 
the WIYN OPTIC imager. Consequently, our survey 
was limited to the galaxy's inner ~2 scale lengths, a 
regime previously probed (in other, more distant galax- 
ies) with integrated - light absorption line spec troscopy 
(e.g. , iBottemal H99l iGerssen et all I1997L I2000D . Since 
these studies have generally found that the mass-to-light 
ratio of inner spiral disks is constant with radius, it is 
not surprising that we find a similar result. 

The exponential curve shown in Figure [¥] is de- 
rived from the B, V and I surface photometry o f 
iButa fc McC^H (fl999h and iWhite fc Botbrml (pol . 
which shows that, over the region of our survey, the disk 
of IC 342 is well fit by a single exponential with a scale 
length of 250" (~4.2 kpc). (There is a possible steepen- 
ing of the profile at R > 700" but that is well beyond our 
frame limit.) Our inferred values of a z follow this curve 
very well, yielding a dynamical scale length (~5 kpc) 
that is not significantly larger than that of the /-band 
light. Moreove r, all our values for a z lie safely under the 
iToomrel ([1964) stability curves of Figure [9l suggesting 
that our estimate of h z is not unreasonable. 

Our most-likely value for the central z-direction ve- 
locity dispersion of IC 342 is a z (0) — 41 ± 6 km s . 
If we assume that the disk T is constant, and adopt 
a scale height of h z = 500 pc, then through equa- 
tion ([1]), the galaxy's central disk surface mass density 
is S(0) = 146 ± 45 {h z /500 pc)" 1 M Q pc" 2 . A compari- 
son of this number with the gala xy's central disk surface 
brightness (/i/(0) = 19.6 ± 0.2: IWhite fc Bothiml 12001 
then yields an /-band disk mass-to-light ratio of Tj ~ 
0.25 ± 0.10 (W500 p c)- 1 for E(B - V) = 0.558 ± 0.069 
(jSchlegel et alill998l ) and a lCardelli etaTI (|1989h redden- 
ing law. Note that our error includes the uncertainties 
in the fit, the surface photometry, and h z (estimated at 
~20%), but not the systematic error associated with the 
disk's thermal structure. As described in §3, this error 
can be as high as ~30%, but it is most likely a systematic 
offset that affects all our measurements. 

8.2. M74 (NGC 628) 

M74 is the most face-on galaxy in our sample, and the 
observed line-of-sight velocity dispersion is almost en- 
tirely due to a z . Like IC 342, multi-band surface photom- 
etry confirms that the disk is well fit b y a single exponen- 
tial over the entire range of our survey (|Natali et al.lll992l : 
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FlG. 9. — Likelihood probabilities for IC 342, as a function of cr s 
(y-axis) and (Jr (a:-axis). The darkness is directly proportional to prob- 
ability, and the white crosses indicate our most likely solutions, without 
any consideration for stability. The solid curves illustrate the maximum 
value of a z , a bove which th e disk becomes unstable to axisymmctric 
perturbations ( Toomrc 1964), and correspond to disk scale heights of 
h z — 400, 500, and 600 pc (from bottom to top). The figure illus- 
trates that our estimates of h z are reasonable, while smaller values are 
problematic, due to stability concerns. 



iKoopmann et al.ll2006[ ) . The galaxy, does, however, have 
a slightly earlier Hubble type (Sc) than IC 342, and is 
slightly smaller: at Paper I's distance of 8.6 Mpc, M74's 
.R-band disk scale length (~3. 2 kpc) is ^25% shorter 
than that of the former galaxy (|Mollenhofi1l2004D . M74's 
maximum rotation velocity is also smaller than that of 
IC 342. 

As seen in Figure a constant h z , constant T model 
is an excellent fit to M74's velocity dispersion. The 
most likely scale length for cr 2 , 3. 3l ; 5 kpc, is sta- 
tistically indistinguishable from that derived from the 
/?-band light, and, with the possible exception of the 
outermost bin, the probabilit ies of Figure [10] are all 
well below the IToomrel ()1964l ) stability curve. In ad- 
dition, the disk's central velocity dispersion, as ex- 
trapolated from the a z profile (53 ± 5 km s" 1 ) is in 
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FlG. 10. — Likelihood probabilities for M74, as a function of a z (y- 
axis) and aji (a:-axis). See the caption of Figurc[9]for more information. 
The three curves show the stability limits for h z — 300, 400, and 500 pc 
(from bottom to top); solutions above these curves are excluded by the 
IToomre] d!964D condition. Within ~10 kpc, a constant scale height of 
~400 pc is reasonable, but for the final bin, a larger value of h z may 
be needed. 

good agreement with the ^60 km s _1 value found from 
integral-held absor ption-line spectrosco py of the system's 
central kiloparsec (jGanda et al.l [20061 ) . Our measure- 
ment, coupled with our adopted scale height estimate 
of h z — 400 pc, implies a central surface mass of 
E(0) = 3051 ^/^/400 pc)' 1 M Q pc" 2 . The surface 
photometry of iMollenhoffl ((2004) combin ed with a fore- 
groun d extinction of E(B — V) — 0.07 (Schlegel ct al. 
119981) then yields T R ~ 1.4 ± 0.3 (h z /400 pc)" 1 , or 
T 7 ~ 1.0 ± 0.3 (fc*/400 pc)- 1 . The fact that this mass- 
to-light ratio is larger than that of IC 342 is consistent 
with the galaxy's earlier Hubble type. 

8.3. M101 (NGC 5^57) 

M101, a grand design Scd spiral, is the largest galaxy 
in our sample, with an i?-band scale length of ^5 kpc 
(|Koopmann et aT| [2b06) , and a maximum rotation speed 
of ^250 km s _1 . Unfortunately, it is also the galaxy 
with the poorest velocity sampling, with only 60 PN 
measurements, including 12 objects whose velocities may 
be compromised by nearby emission. Nevertheless, since 
the system's inclination is only ~17°, the probabilities 
produced by our analysis (Figure ITT]) are relatively flat, 
and the corresponding values of a z are reasonably well- 
defined. Interestingly, the dispersion profile derived from 
our measurements looks nothing like that of a constant 
h z , constant T disk. As Figure [8] illustrates, the best- 
fitting scale length for er 2 is ^5 times larger than that of 
the galaxy's luminosity. 

To explain this result, we first consider the possibility 



that the scale height, h z , is changing over the area of 
our survey. Figure [TT] presents no evidence for this ei- 
ther way, as the probabilities derived from the PN veloc- 
ities are gener ally below the constraint imposed by the 
iToomrel (|1964) condition. However, in order to mimic 
a constant T disk, h z would have to increase by a full 
order of magnitude between ~3.5 kpc and ~18 kpc. 
Such a gradient is unlikely, as it would make M101 
unique: none of the large, ed ge-o n spirals observed by 
Ivan der Kruit &: Searld (|1982f ) and lBizvaev &: Mitronoval 
(|2002ft have such a profile. 

Alternatively, we can assume that h z is constant, 
and ascribe the flat velocity dispersion to a change 
in the galaxy's disk mass-to-light ratio. M101 has a 
central i?-band surfac e brightness of (J,r(0) *~ 20.29 
(jKuchinski et alJ [2000h . and alm ost no foreground ex - 
tinction {E(B - V) = 0.009; ISchlegel et all 119981) . 
These numbers, combined with the most-likely value 
for the disk's central velocity dispersion, cr z (0) — 
38^4 km s _1 , imply a central disk surface mass of S(0) = 
lOOtgJ {h z /600 pc)" 1 M Q pc" 2 , and a central disk mass- 
to-light ratio of T R ~ 0.6±0.3 (h z /600 pc)" 1 . These val- 
ues are reasonable: the surface mass is co nsistent with 
that inferred from the Hi rotation curve ([Walter et al.l 
l2008h . and T falls midway between that of the blue, late- 
type spiral IC 342, and that for the earlier M74 system. 
However, when we adopt these numbers and use the ob- 
served i?-band disk scale length of ~5 kpc, the velocity 
dispersion at ~17 kpc implies a mass-to- light ratio of 
Tr ~ 9 (h z /600 pc)" 1 . This is an extremely large value 
for an actively star-forming disk. 

Again, we note that due to the limited number of plan - 
etary nebulae discovered in the iFeldmeier et alJ ([1997D 
survey, the velocity field of M101 is poorly sampled. 
With only four bins of 15 PNe each, there is not enough 
information to model the galaxy with both a variable h z 
and changing T. Moreover, with such a limited num- 
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FlG. 11. — Likelihood probabilities for M101, as a function of a z (y- 
axis) and <7r (ic-axis). See the caption of Figure[9]for more information. 
The three curves show the stability limits for h z — 500, 600, and 700 pc 
(from bottom to top). As with M74, the data for the inner bins are 
consistent with the nominal scale height, but in the outer disk, a larger 
scale height might be needed. 
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ber of objects, our analysis is susceptible to the effects 
of disk substructure. For example, if some of the PNe 
targeted in our survey actually belong to the remnants 
of a disrupted companion galaxy (like the Milky Way's 
Sagittarius dwarf), then our derived disk mass measure- 
ments would be in error. The only definite statement we 
can make about M101 is that further observations are 
warranted. 

8.4. M94 (NGC 4736) 

The Sab spiral M94 is the earliest galaxy in our sample, 
with a small, but high surface-brightness bulge that in- 
hibits PN detections in the system's central ~0.5 kpc. 
Nevertheless, bec ause the galaxy is re latively small 
(h R = 1.22 kpc; IGutierrez etafl 120091: Ide Blok et all 
l2008f l. nearby (~4.4 Mpc), an d has little foregrou nd ex- 
tinction (E(B - V) = 0.018: ISchlegel et al.l fl998). it is 
the system for which we have the best data (median PN 
velocity uncertainty of ^3 km s _1 ) with a wide spatial 
sampling (out to ~7 disk scale lengths). The result is the 
high-quality dispersion profile seen in Figure which 
clearly shows an exponential decline over the galaxy's 
inner regions, followed by a flattening which begins at 
~4 disk scale lengths. 

Unlike the galaxies of the previous sections, M94 
does not have a simple surface brightness profile: while 
the galaxy's inner disk follows an exponential law, its 
outer regions ar e dominated by the luminosity of a 
faint ste llar ring dSandage 1 119611 : ISchommer fc Sullivan! 
11976ft . iDe Blok et al.l 1)20081 ) have recently analyzed 
M94's profile using 3.6 (urn surface brightness measure- 
ments fr om the Spitzer Infrare d Nearby Galaxy Survey 
(SINGS; iKennicutt et al.l l2003f) . Their data show that 
the galaxy's disk is best modeled by an inner, rapidly 
falling exponential (scale length hjii = 1.22 kpc), which 
is abruptly replaced at R ~ 5 kpc by an almost constant 
surface-brightness (scale length I ir. = 7.16 kpc ) enve - 
lope. The surface photometry of IGutierrez "elaTl (2009) 
confirms this result: in the i?-band, M94 has a rapidly de- 
clining inner disk (/iR,j = 1.22 kpc), a uniform brightness 
outer ring (/i#. = 13.2 kpc between 4.6 and 7.0 kpc), and 
a slowly declining outer disk (h.R. x = 2.98 kpc) . Thus, 
in the classification scheme of lErwin" "eTall (f2005h . M94 
is a Type III "antitruncated" spiral, whose disk surface 
brightness declines with one slope in its inner regions, 
and a different, shallower slope at larger radii. Roughly 
30% of all spirals, and ~50% of all ea rly-type spirals fall 
into this category ([Pohlen et al.l l2008). Furthermore, our 
PN number counts support the classification of M94 as a 
Type III object, as there are far too many planetary neb- 
ulae at large radii to be explained by a single exponential 
law. 

While the cause of the break in Type III disk galax- 
ies is unclear, one possibility is that it is due to 
the i ncreased importance of the thick disk at larg e 
radii l|Yoachim fc Dalcantonl 120061 : iPohlen et aHl2007l ). 
Our kinematic study is consistent with this hypothe- 
sis. Since early-type spiral s typically have h^Jh^ 10 
( Biz yaev fc Mitronovall2002l : Ide Grijslll998l : iKregel et all 
2002), M94's disk scale height should be relatively large, 
h z ~ 300 pc. Indeed, as Figure fl~2l illustrates, a 300 pc 
scale height satisfies the lToomrd ([1964ft stability criterion 
throughout the inner regions of the galaxy, while values 
much smaller than ^200 pc do not. However, outside 



^4 kpc, no thin-disk value of h z is stable; at these radii, 
significantly larger scale heights, even as high as 900 pc, 
are needed. In M94's outer regions, the light we see must 
be coming from a thick disk. 

Based on these data, we decided to model the galaxy 
using two separate kinematic populations. Within 
3.8 kpc, we assume that most of the PNe are part of 
the galaxy's thin disk with h z ~ 300 pc; outside of 
5 kpc, we allow the thick disk to dominate, and set 
the scale height to 3 times that of the thin disk, i.e., 
900 pc. Su ch a ratio agrees with obse rvations of the 
Milkv Wav (jLarsen fc Humphreys! l2003h and NGC 891 
([Morrison et al.lll997ft. and is consistent w i th the range 
of values found by lYoachim fc Dalcantonl (|2006ft for a 
sample of 34 edge-on galaxies. In between these two ex- 
tremes, we define a transition zone, where the effective 
scale height increases linearly with radius. 

The gray curve in the lower panel of Figure [7] illus- 
trates our thin disk-thick disk model, with each com- 
ponent having a fixed mass-to-light ratio. This curve 
is an excellent fit to the data. The central disk ve- 
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FlG. 12. — Likelihood probabilities for M94, as a function of a z (y- 
axis) and <7r (ic-axis). See the caption of Figurc[9]for more information. 
The three curves show the stability limits for h z — 100, 200, and 300 pc 
(from bottom to top). In the inner five bins, a scale height of h z — 
300 pc appears consistent with the data, but at larger radii, a much 
larger scale height is required. 
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locity dispersion implied by our fit is 91 ± 15 km s _1 ; 
when combined with M94's central J? -band disk sur- 
face brightness (/Ar(0) = 18.88 ± 0.05; iGutierrez et al.l 
20091) and foregrou nd extinction (E(B - V) = 0.018; 



Schlegel et al.lfl998l ). this yields a thin-disk mass-to- light 
ratio of T R = 1.8 ± 0.4 {h z /300 pc) -1 . Our thick- 
disk mass-to-light ratio based on this same i?-band sur- 
face photometry and the outer three bins of Figure [7] is 
Tr = 2.0 ± 0.4 (ft z /900 pc) -1 . These values are reason- 
able, given the galaxy's early Hubble type. 

8.5. M83 (NGC 5236) 

The barred spiral M83 is the galaxy for which we have 
the widest spatial coverage (~25 kpc) and the largest set 
of high-precision PN velocities (162). Unfortunately, it 
is also the object with the most puzzling dispersion pro- 
file. As Figure [6] illustrates, except for its innermost bin, 
M83's a z values remain virtually unchanged from ^6 to 
^20 kpc, a region which spans ^4 disk scale lengths. 
One reason for this behavior may be the presence of the 
system's strong bar, as the structure's outer Lindblad 
resonance is located between ou r second and third ve - 
locity bins (at 290" or ~6.83 kpc; iLundgren et al.ll200l . 
Since all the stellar velocities in the area are affected, it 
is quite possible that our derived value for the region's 
a z is an underestimate. These points should therefore be 
excluded from the analysis. 

The rest of M83's dispersion profile is not as easy 
to interpret. Unlike M94, M83 shows no evidence for 
a Type III anti-truncated disk or a flattened surface- 
brightness profile at large radii. Indeed, according to 
de Jong (2008, private communication), the galaxy's H- 
band surface brightness profile declines with a single ex- 
ponential scale length of 80" (1.88 kpc) over the entire 
region of our survey. Thus, the relatively large velocity 
dispersion in the outer disk cannot be explained by an 
excess of stellar mass in the region. 

Another puzzle comes from the galaxy's rather small 
radial scale length. Although M83 is the fastest rota- 
tor in our sample (with a v max of 255 km s -1 ), it also 
has the second smallest value of Iir: 2.45 kpc in R and 
1.88 kpc in H (iKuchinski et al.ll2000l de Jong 2008, pri- 
vate communication). If the scale length to scale height 
ratio is similar to that of other galaxies, i.e . , hp /h z ~ 10 
dBizvaev fc Mitronova|[2002t Ide Grilill998h iKregei 1 et all 
l2002h . then this implies a value for h z near ~200 pc. 
However, as Figure [13] demonstrates, scale heights as 
small as this yield solutions for a z that ar e unstable to 
axisymmetric perturbations (|Toomre!ll964h . As a result, 
any attempt to fit the data by varying h z while keeping 
hji constant leads to either an unphysically large scale 
height, an exceptionally low value of hn/h z , or a large 
central velocity dispersion that is completely inconsistent 
with the galaxy's rotation curve. 

Given these difficulties, we chose to model M83 with 
a two component, thin disk-thick disk model. For the 
galaxy's thin disk, we adopted a scale length of 4.0 kpc, 
a central velocity dispersion of <r 2 (0) ~ 73 km s , and 
a scale height of 400 pc; for the systems' thick disk, we 
assumed Hr ~ 10 kpc, a z (0) — 40 km s _1 , and h z = 
1200 pc. As the curve in the lower panel of Figure [6] 
indicates, this type of model, although arbitrary, does 
do a good job of reproducing the observed behavior of 
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FlG. 13. — Likelihood probabilities for M83, as a function of er s (y- 
axis) and itr (a;-axis). See the caption of Figure[9]for more information. 
The four curves show the stability limits for h z — 200, 300, 400, and 
500 pc (from bottom to top). As in M94, the data suggest that larger 
scale heights arc needed to maintain stability against axisymmetric 
perturbations (Toomrc 1964). 

If we adopt this model, and use a central f?-band 
disk surface brightness of /J.(R) = 19.07 (jKuchinsk i et al.l 
2000) and a f oreground Galactic e xtinction of E(B — 
V) = 0.066 (|Schlegel et all fl998l) . then M83's cen- 
tral disk mass-to-light ratio becomes Tr = 1.0 ± 
0.5 (/i z /400 pc) -1 , where the 50% uncertainty is a con- 
servative estimate. In keeping with the galaxy's Hubble 
type, this value is intermediate between the lower values 
of IC 342 and M101 and the high values for M74 and 
M94. Of course, because the galaxy's surface brightness 
declines with radius while its z velocity dispersion does 
not, the inferred value for T increases dramatically in 
the galaxy's outer regions. 
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9. DISCUSSION 

Except for M101, the data from Figures[3]through[5]of- 
fer a consistent picture for the mass profile of spiral disks. 
Within the inner ~3 to 4 disk scale lengths, the ratio of 
o%/ti(R) in IC 342, M74, and M94 is flat with radius. 
Through equation ([1]), this implies that the stellar ver- 
tical velocity dispersion is consistent with that expected 
from a constant h z , constant T disk. (The data for M83 
can also be interpreted in this way, if one accepts that 
the stellar kinematics between ~6 and ~7 kpc are af- 
fected by the bar's outer Lindblad resonance). In the 
latest-type spirals, IC 342 and M101, the disk mass-to- 
light ratio is less than one; in the earlier systems, Tr is 
roughly between one and two. Although these numbers 
have a degree of systematic uncertainty due to our lack 
of knowledge about the disks' vertical structure, they 
are nonetheless in the range expected for normal, star- 
forming stellar populations. 

Our measurements of the disk velocity ellipsoids are 
quite uncertain. Since the galaxies in our sample were 
selected on the basis of their low inclination, our ability 
to measure aR, the disk velocity dispersion in the radial 
direction, is very limited. Consequently, as Figure 1141 
illustrates, a marginalization over the likelihoods of Fig- 
ures |9][T3] places few constraints on a z /aR, especially in 
our two most face-on galaxies, M74 and M101. The trend 
toward higher a z /aR in the three more inclined systems 
is likely a numerical artifact caused by switching from 
&R-&Z space to (<J z /an)-a z space for this analysis. 

Finally, the data in Figure Q3] can be used to confirm 
our estimates of each galaxy's asymmetric drift. The full 
expression for the drift velocity in a stellar disk is 
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where v is the stellar density (|Binnev fc Trema inc 1987). 
However, if the disk is roughly isothermal with constant 
h z , then the epicyclic approximation gives 

V 2 - {V*) 2 ~2V c v asd (13) 
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Finally, Figure [14] suggests that the shape of a disk's 
velocity ellipsoid does not change much with radius. In 
this case, the application of equation |T]) creates a very 
simple expression for asymmetric drift 
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where we have substituted in Hr by assuming a constant 
mass-to-light ratio in the stellar disk. 

This equation yields roughly constant asymmetric drift 
velocities throughout the three disks where our assump- 
tions appear valid. In the inner regions of our earliest spi- 
ral, M94, equation (fl"4"|) implies an asymmetric drift ve- 
locity of v as d ~ 40 km s _1 , exactly the value assumed in 
§4.2. Conversely, in the two latest galaxies in our sample, 
our derived values for the drift velocity, v as d ~ H km s _1 
(for IC 342) and -8 km s _1 (for M101) are slightly lower 
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FlG. 14. — Relative likelihoods for <7 z /i7r in the five disks observed in 
this survey. The PNe in each galaxy were divided as equally as possible 
into three regions (inner, middle, and outer). Because an is poorly 
constrained in these low-inclination systems, the likelihood curves are 
all very broad. In our moderately inclined galaxies (IC 342, MS3, and 
M94), the results seem to favor higher values of (T z /cth, but this could 
be an artifact caused by switching from (7r-ct z space to {a z / <tr)-(t z 
space. There is no evidence for a shape gradient in any of the galaxies. 



than the values assumed above (20 and 25 km s _1 , re- 
spectively). Of course, given the uncertainties associated 
with our estimates of aR, the numbers are still consistent. 
Finally, the face-on nature of M74 and disturbed rotation 
curve of M83 (due to its bar), preclude statements about 
the asymmetric drifts in these galaxies. However, it does 
appear that our assumption of asymmetric drift being 
between 10% and 20% of the rotation speed is valid. 

9.1. Inner Disks 

Most analyses of the rotation curves of spirals galaxies 
make use of the maximal-disk hypothesis, wherein the 
rotation curve near the center of the galaxy is assumed 
to be defined entirely by the potential of a thin disk. Our 
measurements of a z allow us to test this assumption and 
probe the distribution of the systems' dark matter halos. 
The circular velocity produced by an exponential disk 
with radial scale length, h R , and central surface mass, 
S , is given by 

V C 2 (R) = 4nGXoh R y 2 [I (y)K (y) - (15) 

where y = R/{2hn) and I n and K n are modified Besse l 
functions (jFreema 3 [19701 IBinnev fc Tr emaind Il987h . 
When combined with equation |T]), this implies a relation 
between the disk's rotation speed, ^-velocity dispersion, 
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FlG. 15. — The rotation curves for the inner disks of our fi ve program gal a xies. The solid black curves display the general motion of the Hi gas 
(the black dots in M94 show the THINGS rotation curve from lde Blok et al.l J2008D 1: the green curves show the contribution of the disk under the 
assumption of an intermediate {scc\i{z f h z )) vertical structure and the values of h z given in the last column of Table[2] The blue curves show the 
contribution of the residual halo component, with the NFW profiles shown with long orange dashes and the pseudo-isothermal models denoted by 
short red dashes. The bulge of M94 (indicated by a dotted dark red curve) is a substantial and uncertain contributor to the mass, so no residual 
halo was fit to this system. (See the on-line article for a color version of this figure.) 



and scale length to scale height ratio, i.e., 

v c 2 (R) 
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where K — 1.70517T for a disk with a vertical distribution 
intermediate between exponential and isothermal. Since 
the rotational velocities o f our disks are we l l measured 
(ICrosthwaite et al.l 12000' iFathi et al.l 120071 : I Park et all 



120011: Ide Blok et all 120081: IZaritskv et al.l fl990f i . and 

the behavior of hn/h z with Hubble type a nd rota- 
tion s peed is know n fro m the literature (e.g., Ide Griisl 
119981: iKreeel et all l2002t lYoachim fc Dalcantonl I2006D . 
our measurements of a z (0) place an immediate constraint 
on the importance of dark matter near the centers of 
galaxies. 

The result of this analysis is illustrated in Figure [T5] 
and summarized in column 7 of Table [2j In brief, the five 
disks examined in this survey display a range of proper- 
ties, depending on the Hubble type of the system. The 
disks of our latest galaxies, IC 342 and M101, are clearly 
sub-maximal, requiring scale-length to scale-height ra- 
tios hft/h z > 40 to attain the maximal condition. Un- 
less these disks are extremely thin (100 pc for IC 342 
and 60 pc for M101), their dark matter halos must be 
contributing significantly to the potential of their inner 
regions. In the Sc galaxies M74 and M83, the situa- 
tion is not as clear. For M74's disk to be maximal, 
hR/h z ~ 18 (i.e., h z ~ 170 pc). Although this seems 
significantl y larger than t he value of ~8 expected for an 
Sc galaxy (|de Griislll998D . the two values can be forced 
into marginal agreement if one assumes that the vertical 
structure of the disk is exponential, rather than interme- 
diate or isothermal. Similarly, M83 has a requisite scale 
ratio near the edge of the range expected for SBc galax- 



ies, h,R/h z ~ 10. Only the Sab galaxy M94, with its 
high surface-brightness bulge has a disk that can be con- 
sidered maximal: if we model the galaxy's in ner regions 
with a n h z ~ 300 pc disk and a T# ~ 1.3 Hernquist 
(1990) bulge (a — 0.586 kpc), then the rotation curve 
of the central ~1.5 kpc is entirely explained by baryons. 
We note that with this simple model, the mass-to-light 
ratio of M94's inner disk is actually greater than that 
of its bulge; however, if the disk's vertical structure is 
modeled as isothermal, rather than intermediate, or if 
the disk's scale height is smaller than ^300 pc, then the 
mass-to-light ratios come into much better agreement. 

Since four of our galactic disks are sub-maximal, it is 
a simple matter to remove the contribution of their rota- 
tion curves and examine the potential of the remaining 
dark matter component. When describing halo distribu- 
tions, two models are generally considered: one following 
the predictions of iNavarro et al.l (|1996t 119971 hereafter 
NFW) and one obeying a pseudo-i sothermal (ISO) law 
(jBahcall fc Soneirdl 980 iKentl 198(1 . The rotation curve 
for the former is given by 



Vnfw(R) — Vi 



200 t 



/ln(l + cx) — cx/(l + cx) 
x[\n{\ + c) - c /(l + c)] ' 



(17) 



where x = i?/i?200j -R200 is the radius at which the den- 
sity contrast over the background exceeds 200, V200 is the 
circular velocity at i?2oo > and c is a concentration param- 
eter, which physically should always be greater than 1. 
In contrast, the rotation curve for the pseudo- isothermal 
law is 



Viso(R) = JinGpoRl 



1 R c t ( R 
1 arc tan 

R \Rc 



(18) 
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where po is the central density of the halo and Rq is 
the core radius. Figure [TBI fits both these profiles to our 
residual rotation curves. As the figure illustrates, the 
isothermal model fits the data extremely well with mean 
variances of 0.24, 0.17, 0.04, and 4.52 km 2 s" 2 for IC 342, 
M74, M101 and M83, respectively. (In M83, the lack of 
PN identifications in the central ~1 disk scale length lim- 
its the accuracy of our extrapolation.) In contrast, the 
NFW fits are all very poor, with large mean variances 
(39.2, 37.2, 8.51, and 174.8 km 2 s~ 2 for these same four 
galaxies), unphysical concentration indices (c <C 1) and 
density distributions that overpredict the inner rotation 
curve. As in the case of low-surfac e brightness galax- 
ies fe .g..lKuzio de Narav et alll2008l ) and dwarf galaxies 
fe.g.. lOh et alJl2008ft . the inner rotation curves of normal 
spirals are not well represented by NFW profiles. 

9.2. Outer Disks 

In contrast to the galaxies' inner regions, the outer 
disks of spirals are not constant T, constant h z systems. 
As the data of M83, M94, and perhaps M74 show, in 
these distant regions, it is the product of surface mass 
times scale height that appears constant. Specifically, in 
M94, the ratio of <t 2 /(G/j) = KTh z has one value in the 
galaxy's inner regions, and another larger value in the 
outer disk, and in M83, the ratio blows up at large radii. 
The same effect can b e seen in the outer regions of M33 
Ciardullo et al.|[2004f ). though the larger inclination of 
this galaxy creates an ambiguity in the interpretation.) 
At radii greater than ^3 disk scale lengths, the disks of 
spiral galaxies arc cither flaring, increasing their mass- 
to-light ratio, or changing their luminosity profile from 
that of a single scale-length exponential. 

To some extent, all three of these phenomena must 
be happening. Since many galaxies have large Hi en- 
velopes, an increase in the disk mass-to-optical light ratio 
at large radii is not unexpec ted. However, meas urements 
from the THINGS survey ([Walter et all 12008( 1 demon- 
strate that this effect is minor: even in the outer re- 
gions of M83 and M94, Hi contributes no more than 
~1Q% of the kinematically-estimated disk mass. Simi- 
larly, many spirals do not have simple exponential disks 
in their outer regions. In M94, a shallo wing of the sur- 
face b rightness pro file was recorded by iGutierrez et al.1 
([2009T1 and SINGS (|de Blok et al.l[2008h . and this excess 
luminosity goes a long way towards explaining the sys- 
tem's mass-to-light ratio. Unfortunately, not all spirals 
have such breaks, and the disk of M83 appears to follow 
the same exponential over the entire range of our survey. 
Finally, stability arguments demand that some flaring 
occurs in the outer parts of spirals. In all of our galaxies 
(except IC 342, whose data only extends to ~2 disk scale 
lengths), h z likely incr eases past ~3 d isk scale lengths, 
in order to satisfy the iToomrei (jl964f ) criterion against 
axisymmetric perturbations. 

Is there any direct evidence for flaring in the outer 
disks of spirals? Studies of edge-on galaxies have gener- 
ally found that galactic scale heights are roughly con- 
stant with radius (e.g.. | v an der Kruit fc Searld fl98l 
iBizvaev fc Mitronoval [2002). Even when radial gradi- 
ents have been found, they have generally been small, 
and linked to the presence of the thic k disks of early- 
type systems ([de Griis fc Peletier1ll997t l. However, these 
studies have mostly been confined to the systems' bright 



inner regions. In NGC 891, which has extremely deep 
surface photometry derived from individual star counts, 
a thick, highly-flattened, low surface brightness structure 
does extend well bey ond the limits of the galaxy's thin 
disk ([Mouhcindl2009[ l. 

Mechanisms do exist which can cause the scale height 
of star s to increase rapidly at large radii. iNaravan fc Jod 
(2002) have shown that self-gravitating disks acting un- 
der the influence of Hi and H2 gas and a dark matter 
halo can develop moderately large scale height gradients 
in their outer regions (R > 2 disk scale lengths). Al- 
ternatively, flared-disks and higher than expected val- 
ues of (j 2 can be caused by in teractions with halo sub - 
structures. In fact, Model F bv lHavashi fc Chibal (|2006(1 . 
which simulates the heating of the Milky Way's disk by 
a population of ^300 moderately massive (~10 8 ' 5 M Q ) 
dark matter subhalos, produces a velocity dispersion pro- 
file t hat is an excellent matc h to that seen in M83 and 
M94 (H errmann et al.|[2009al ). Recent simulations by S. 
Kazantzidis (2009, private communications) confirm this 
result, as they reproduce the full <j z profiles of the two 
galaxies extremely well. Thus, the flaring inferred by our 
kinematic data is certainly plausible. 

10. CONCLUSIONS 

An analysis of 550 precise PN velocities has yielded 
kinematic mass estimates throughout the disks of five 
low-inclination spiral galaxies. In the inner regions of 
IC 342, M74, and M94, we find that the velocity dis- 
persion perpendicular to the galactic disk exponentially 
decays with a scale length twice that of the optical light. 
While our observations of the remaining two galaxies 
yield ambiguous results (due to a Lindbland resonance 
in M83 and a poorly sampled velocity field in M101), 
the data are consistent with expectations from a con- 
stant scale height, constant mass-to-light ratio disk. In 
the later-type (Scd) systems, these disks are clearly sub- 
maximal, with surface mass densities less than a quarter 
that needed to reproduce the central rotation curves. In 
earlier (Sc) galaxies, this discrepancy is smaller, but still 
present; only the very-early Sab system M94 has evidence 
for a maximal disk. 

When we subtract off the disk contribution from the 
Hi rotation curve, we infer the potential of the galax- 
ies' invisible dark halos in all our galaxies except M94. 
(The uncertainties in the contribution of M94's bright 
bulge preclude such an analysis in this earliest galaxy 
in our sample.) In all four cases, these halos are best 
fit with a simple pseudo-isothermal model; NFW profiles 
are poor fits to the data, as they always overpredict the 
residual rotation curve. As in the analysis of low surface- 
brightness and dwarf galaxies, we find no evidence for 
central dark matter cusps in our systems. 

Interestingly, once outside of ^3 disk scale lengths, 
the z-direction velocity dispersion of galactic disks no 
longer declines, but instead stays constant with radius. 
This behavior, coupled with stability arguments, sug- 
gests that the disks of spiral galaxies flare substantially 
at extremely large radii, and may have strongly increas- 
ing mass-to-light ratios. Whether this flaring is related 
to the increased importance of the thick disk, or heating 
of the thin disk by halo substructure, is unclear at the 
present time. To further investigate this phenomenon, 
deep surface photometry of the extreme outer disks of 
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edge-on galaxies is needed. 

Our kinematic mass estimates for galactic disks are 
consistent with expectations. The galaxies with the lat- 
est Hubble types (IC 342, M101) have the smallest mass- 
to- light ratios, while earlier systems (M74 and M94) have 
values of T greater than one. Unfortunately at this time, 
no deep, multi-band optical and IR surface photometry 
exists for any of our galaxies. Once these data are ac- 
quired, we will be able to compare our kinematic values 
for T to stellar mass-to-light ratio estimates from pop- 
ulation synthesis models. Such a confrontation will not 
only produce a better constraint on galactic populations, 
but also will improve our knowledge of the phase-space 
distribution of galactic disks. 
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